AD^ftllO  455  WISCONSIN  UNIV-MADISON  MATHEMATICS  RESEARCH  CENTER  F/G  7 

OYNAMIC  PHASE  TRANSITIONS  IN  A  VAN  DER  WAALS  FLUID. <U> 

NOV  81  M  SLEMROD  DAAG29-80-C-0041 

UNCLASSIFIED  MRC-TSR-2298  NL 


n»«C0M.  «  __  »«I0455 


athematics  Research  Center 
liversity  of  Wisconsin-Madison 
0  Walnut  Street 
adison,  Wisconsin  53706 

November  1981 


Received  September  18,  1981) 


tonsored  by 


Approved  for  public  release 
Distribution  unlimited 


A?-*-  j  oo 


O.  S.  Army  Research  Office 

P.  O.  Box  12211 
Research  Triangle  Park 
North  Carolina  27709 


National  Science  Foundation  Air  Force  Office  of 
Washington,  DC  20550  Scientific  Research 

Washington,  DC  20332 


82  02  03  086 


UNIVERSITY  OF  WISCONSIN  -  MADISON 
MATHEMATICS  RESEARCH  CENTER 


DYNAMIC  PHASE  TRANSITIONS  IN  A  VAN  DER  WAALS  FLUID 

* 

M.  Slemrod 

Technical  Summary  Report  #2298 
November  1981 

/■  ABSTRACT 

—  JP'This  paper  proves  the  existence  of  traveling  wave  solutions  connecting 
liquid  and  vapor  phases  in  a  van  der  Waals  fluid.  The  main  constitutive 
assumptions  are  that  the  fluid  be  an  elastic  fluid  (with  pressure  given  by  the 
van  der  Waals  equation  of  state)  possessing  a  higher  order  correction  given  by 
Korteweg’s  theory  of  capillarity  and  the  fluid  is  a  conductor  of  heat  with 
large  specific  heat  at  constant  volume.  The  main  mathematical  tool  in  the 

analysis  is  the  Conley- Easton  theory  of  isolating  blocks.  ■ 
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SIGNIFICANCE  AND  EXPLANATION 


In  the  study  of  phase  transitions  in  fluids  a  typical  model  given  in 
classical  thermodynamics  is  provided  by  the  equilibrium  configurations  of  a 
van  der  Waals  fluid.  Such  equilibria  show  a  fluid  may  exist  in  two  phases, 
liquid  and  vapor.  Relatively  untouched  in  classical  thermodynamics  is  the 
non-equilibrium  case.  This  paper  attempts  to  study  this  problem  from  the 
viewpoint  of  examining  the  dynamics  of  liquid  -  vapor  phase  transitions.  The 
main  mathematical  idea  is  to  look  for  traveling  wave  solutions  of  the  relevant 
balance  laws  and  show  the  desired  connections  from  one  phase  to  another  can  be 
made.  Since  the  equations  are  analogous  to  those  modeling  phase  transitions 
of  the  "martensitic"  or  "shape  memory"  type  in  solids  the  results  given  here 
may  prove  of  value  in  that  case  as  well. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


DYNAMIC  PHASE  TRANSITIONS  IN  A  VAN  DER  WAALS  FLUID 


M.  Slemrod* 

0.  Introduction 

The  purpose  of  this  paper  is  to  prove  the  existence  of  traveling  wave 
solutions  connecting  liquid  and  vapor  phases  in  a  van  der  Waals  fluid.  In  an 
earlier  paper  [1]  I  considered  this  problem  when  the  fluid  is  emersed  in  a 
heat  bath  thus  forcing  temperature  to  be  constant.  Here  no  such  assumption  is 
made;  the  motions  are  not  isothermal.  The  constitutive  hypotheses  made  are 
that  the  fluid  be  (i)  an  elastic  fluid  with  pressure  given  by  the  van  der 
Waals  equations  of  state  (see  (1.7))  possessing  (ii)  a  higher  order 
correction  term  given  by  Korteweg's  theory  of  capillarity  (see  (1.5))  and  that 
the  fluid  be  (iii)  a  heat  conductor  with  a  large  specific  heat  at  constant 
volume.  Under  these  hypotheses  I  show  that  given  an  equilibrium  state  in  the 
stable  part  of  the  liquid  (respectively,  vapor)  phase  there  is  a  second 
equilibrium  in  the  vapor  (liquid)  phase  to  which  it  may  be  connected  with  a 
traveling  wave  of  positive  (negative)  speed.  Thus  all  the  vapor  (liquid)  is 
converted  into  liquid  (vapor).  I  also  prove  that  in  both  cases  the  liquid 
equilibrium  phase  is  at  a  higher  temperature  then  the  vapor  equilibrium  phase 
to  which  it  is  connected. 

The  main  tool  of  the  analysis  is  the  Conley-Easton  theory  of  isolating 
blocks  [2].  This  theory  has  been  applied  by  Carpenter  [3]  to  prove  the 
existence  of  traveling  wave  solutions  to  a  generalization  of  the  Fitzhugh- 
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Nagumo  equations  modeling  nerve  impulse  transmission.  Surprisingly  the 
equations  governing  traveling  wave  solutions  in  the  phase  transition  problem 
given  here  and  the  generalized  Fitzhugh-Nagumo  equations  are  similar.  In  fact 
I  show  a  modification  of  Carpenter's  approach  to  the  wave  impulse  equations 
yields  the  desired  traveling  wave  solution  to  the  phase  transition  problem. 

While  this  paper  is  devoted  to  fluid  mechanics,  it  is  closely  related  to 
a  recent  paper  of  James  [4].  In  that  paper  he  used  a  similar  constitutive 
model  to  discuss  phase  transitions  in  a  solid  exhibiting  shape  memory.  He 
showed  that  a  knowledge  of  the  non-isothermal  dynamical  phase  transitions  can 
be  used  to  determine  certain  constitutive  parameters,  e.g.  the  well  known 
Maxwell  line.  Since  the  theory  given  here  gives  conditions  for  such  dynamical 
phase  transitions,  it  may  prove  useful  in  the  study  of  shape  memory  solids  as 


well  as  fluids 


1.  One  dimensional  Lagrangian  description  of  compressible  fluid  flow. 

We  follow  the  presentation  of  Courant  and  Friedrichs  [5]  of  a  Lagrangian 
description  of  compressible  fluid  flow  based  on  the  law  of  conservation  of 
mass.  The  fluid  flow  is  thought  of  as  talcing  place  in  a  tube  of  unit  cross 
section  along  the  x-axis.  We  attach  the  value  X  »  0  to  any  definite  "zero" 
section  moving  with  the  fluid.  For  any  other  section  we  let  X  be  equal  in 
magnitude  to  the  mass  of  the  fluid  in  the  tube  of  unit  cross  sectional  area 
between  that  section  and  the  zero  section.  Analytically  the  quantity  X 
satisfies  the  relation 

X  *  ^x(O^t)  p(x't)dx  •  (1.1) 

Here  p(x,t)  denotes  the  density  at  position  x  and  time  t  and  x(X,t) 

denotes  the  position  of  a  particle  which  encloses  a  mass  X  of  fluid  in  the 
tube  bounded  by  x(X,t)  and  x(0,t).  Differentiation  of  (1.1)  implies  1  * 
xx<X,t)p(x(X*t),t).  Set  p(x(X,t),t)  -  P(X,t),  w(x,t)  »  p(X,t)-1  (the 
specific  volume),  u(x,t)  ■  xfc(X,t)  (the  velocity). 

Also  we  let 


the  pressure, 
the  stress, 

specific  internal  energy, 
specific  total  energy, 
specific  heat  absorption, 
heat  flux, 

specific  body  force, 

absolute  temperature, 

specific  entropy, 

specific  Helmholtz  free  energy, 

specific  heat  at  constant  volume, 

coefficient  of  thermal  conductivity. 
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Specifically  this  means 


T  -  -p(w, 0 )  +  M(w,9)ux  +  D(w,0)wx  -  c<w*6)wxx  » 
e  =  e(w,0)  , 

h  =  h(w,ux,0,9x)  ,  (1.5) 

Hi  =  4»(w,0)  , 

n  =  n(wf0)  , 

where  n,p  satisfy  the  thermodynamic  relations 

n  *  "  P  -  '  '  (1*6) 

and  p,  D,  C  are  small. 

The  reason  for  introducing  Korteweg’s  theory  is  the  following:  Our  main 
interest  is  in  the  case  where  p  is  a  non-monotone  function  of  w  for  fixed 
9  so  that  pw  can  be  both  positive  and  negative  for  different  values  of 
(w,0).  For  the  case  of  an  elastic  fluid  t  =  -p(w,0),  (1.4)  is  a  mixed 

hyperbolic-«lliptic  initial  value  problem  on  fixed  isotherms  0  5  constant. 

But  such  a  problem  is  ill-posed  for  initial  data  in  lying  the  in  the  elliptic 
domain.  One  way  to  circumvent  this  difficulty  is  to  consider  higher  order 
effects  which  become  important  precisely  in  the  elliptic  domain,  namely 
viscosity  and  interfacial  capillarity.  One  theory  of  capillarity  is 
Korteweg’s  and  a  discussion  of  his  work  may  be  found  in  the  monograph  of 
Truesdell  and  Noll  [7).  We  note  that  Korteweg's  theory  has  most  recently  been 
reconsidered  by  Serrin  (8]  who  has  applied  it  to  the  study  of  finding 
conditions  for  equilibrium  of  fluid  and  vapor  phases  in  a  van  der  Waals 
fluid.  (A  related  theory  of  capillarity  based  on  statistical  mechanics  has 
been  presented  in  [91.) 

We  shall  assume  p  in  the  constitutive  relation  (1.5;a)  is  given  by  the 
van  der  Waals  equation  of  state 


constants  (10] ,  (11] ,  [12] 


Isotherms  of  p  are  represented  in  Figure  1. 


PA 


b 


Figure  1 

Also  we  note  from  (1.6)  that  $  is  given  by 

♦  -  -R8  in(w-b)  -  ~  +  F(0) 
w 

A  A  A 

where  F(0)  is  an  arbitrary  function  0.  Since  e  ■  i|>  +  0n  it  follows  that 

e(w, 0 )  «  -  £  +  F(0)  -  OF*  t Q )  • 

As  a  simplifying  assumption  we  set  (see  [12],  p.  74)  F(0)  *  -cv84n0  + 

constant  where  cy  is  a  positive  constant,  to  obtain 


2.  Dynamic  phase  transitions 

We  fix  our  attention  on  a  typical  non-monotone  isotherm  as  shown  in 

Figure  1.  The  domains  (b,w  )  and  (wa,*)  will  be  called  the  a-phase  and 

a  p 

6-phase  respectively; 


(i) 

Pw(w,0)  <  0  ,  0  < 

b  <  w  <  wo,  wp  <  w  , 

(ii) 

pw(%V0)  "  pw<W6'0> 

-  0  , 

(iii ) 

p^(w,0)  >0  if  wa 

<  «  <  w&  . 

(Of  course  w  ,  wQ  depend  on  the  choice  of  non-monotone  isotherm. )  The  o- 
a  p 

phase  corresponds  to  the  fluid  being  liquid/  the  6-phase  corresponds  to  the 
fluid  being  vapor. 

Assume  w_  is  a  constant  value  of  the  specific  volume  lying  in  a-phase 
of  the  6  isotherm  as  shown  in  Figure  2.  Our  goal  will  be  to  see  how  a 
homogeneous  equilibrium  state  of  (1.4)  (with  b  =  q  =  0)  w  =  w_,  u  =  u_, 

0  »  8_  associated  with  the  fluid  being  liquid  can  be  dynamically  transferred 
to  (or  from)  a  second  equilibrium  state  w  *  w+,  u  =  u+,  0  *  0+  associated 
with  the  fluid  being  vapor.  Mathematically  this  means  we  shall  try  to  find  a 
traveling  wave  solution  of  (1.4)  connecting  equilibrium  states  (w  ,u  ,0  ) 
and  (w+,u+,0  +  ),  where  w_  is  in  the  a-phase  of  0  isotherm  and  w+  is 
in  the  6-phase  of  the  0+  isotherm. 

In  order  to  simplify  our  calculations  we  take 

M(w, 0 )  ®  v  >  0,  U  a  small  constant  (the  viscosity)  ; 

2 

C(w,0)  -  M  A  ,  A  a  positive  constant  ; 

h  ■  u0^  ,  a  positive  constant  ;  (2.1) 

D  =  0  ; 

b  =  q  =  0  . 
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As  mentioned  above,  we  seek  a  traveling  wave  solution  for  (2.2)  of  the 


x-ut 

form  u  =  u(C),  w  =  w(C),  0  =  0(5),  where  5  =  — - —  (U  the  speed  of  the 
traveling  wave),  subject  to 

w(  -*)  =  w_,  u(-«)  =  u_,  0(-®)  =  0_  .  (2.3) 

Substitution  into  (2.2)  yields  the  system 

« 

-  Uu'  =  ( -p  +  u*  -  Aw") 

—  Uw 1  sb  u 1 

(2.4) 

-  UE'  =  (u(-p  +  u'  -  AW")'  6" 


where 


d 

d5* 


Set  e  * 
-»  to  t.  We 


e(w_,0  ),  p  =  p(w  ,0  ),  E_  *  —  +  e  and  integrate 
then  find  u,  w,  0  satisfy  the  system 


(2.4)  from 


-  U(u-u  )  =  -(p-p  )  +  u'  -  Aw" 


-  U(w-w  )  =  u-u 


(2.5) 


-U(E-E  )  =  u(-p+u'  -  Aw")  +  u  p  +  -2.  0' 


Finally  we  use  (2.4;b)  and  (2.5;b)  to  eliminate  u  and  u'  from  (2.5;  a,c). 
A  straightforward  computation  shows  w, 0  satisfy  the  equations 
w*  =  v 

2 

Av*  =  -Uv  -  U  (w-w_)  -  p(w, 9 )  +  p  (  2.  6; U, y ) 

2 

01  “  a  u{-(e(w,9)  -  e_)  -  p_(w-w_)  +  y-  (w-w_)2}  . 

In  terms  of  (2.6;U,y)  our  goal  outlined  at  the  beginning  of  this 
section  becomes  the  following:  Find  a  solution  of  (2.6;U,M)  connecting  the 
equilibrium  points  w  =  w_,  v  =  0,  0=0,  w  =  w+,  v  *  0,  9  =  9+,  of 
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(2.6jU,m)  where  w_  is  in  the  a-phase  of  the  0_  isotherm  and  w+  is  in 
the  0-phase  of  the  0+  isotherm. 

We  note  (2.6>U,u)  is  remarkably  similar  to  the  system 

W*  -  V 

AV'  >  -UV  +  G(W,0)  (FN  -  TW) 

0*-  uir’fKw.e) 

governing  traveling  wave  solutions  of  a  generalization  of  the  Fitzhugh-Nagumo 
equations  modeling  nerve  impulse  transmission.  One  obvious  difference, 
however,  between  (2.6;(J,ii)  and  (FN  -  TW)  is  that  the  equilibrium  points  of 
(2.6;U,m)  depend  on  U,  in  (FN  -  TW)  they  don't.  Nevertheless  the  two 
sets  of  equations  are  sufficiently  alike  so  that  a  modification  of  Carpenter's 
argument  [3J  proving  the  existence  of  a  heteroclinic  solution  of  (FN  -  TW) 
will  yield  existence  of  a  solution  to  (2.6;U,lO  connecting  a-  and  0- 
phases.  It  is  this  topic  which  will  be  pursued  in  the  next  sections. 
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3.  The  singular  solution 


Since  u  is  assumed  to  be  small,  a  useful  first  step  in  the  analysis  of 
(2.6;U,m)  is  the  study  of  (2.6;U,0).  This  has  been  done  in  [1]  and  we 
recall  those  results  here.  First,  however,  we  introduce  some  notation. 

For  w_,U  (respectively  w+,U)  given  in  the  a-phase  (B-phase)  of  a 
non-monotone  8_  (8  +  )  isotherm  let  w+(U,0  )  (w  (U,0+))  denote  the  solution, 
if  it  exists,  of 

-  p(w,0  )  +  p(w  ,0  ) 


w-w 


(u2  _  -p(vV  +  P(W'V) 


w  -w 
+ 


lying  in  0-phase  (a-phase).  Also  set 

f_(C;U)  =  U2(5-w_)  +  p(t,0J  -  p(w  ,0  )  , 
f .(CjU)  =  U2(c-w  )  +  p(C,0  )  -  p(w  .0  )  . 

We  then  recall  from  [1]  the  following  lemma. 

Lemma  3. 1  ( i )  Let  w_  be  given  in  the  a-phase  of  the  8  non-monotone 

isotherm.  If 

w  (0,0  ) 

(I)  /w  "  fJC;0)dC  <  0 

then  there  exists  a  unique  U*(0_),  0  <  U*(0  )  <  TJ(0  )  so  that  (2.6; 

*  * 

U  (0_),O)  possesses  a  solution  (w,8),  w(-*)  -  w  ,  w(+*»)  -  w+(U  (0  ),0  ), 

0(0  =  0_.  Here  U(0  )  is  such  that 

w+(U(8  ), 0  ) 

L  f  <C>  U( 9  ))dC  *  0  . 


w+(U  lies  in  the  0-phase  of  the  8  isotherm. 


(ii)  Let  w+  be  given  in  the  3-phase  of  the  8+  non-monotone 


isotherm.  If 


<«>  C(0.9  ) 


then  there  exists  a  unique  U  (0+),  U(8+)  <  U+(0+)  <  0,  so  that  (2.6; 

*  * 

U  (0+),O)  possesses  a  solution  (w,0),  w(-“)  =  w  (U  (0+),9+),  w(+»)  =  w^, 

0(0  =  8  .  Here  U(8  +  )  is  such  that 


w 

/  +  „  f  K;  U(  8  )  )dC  -  0  . 

w_(U  (9  +  ),0+) 


* 

w  (U  (®+).®+)  lies  in  a-phase  of  the  8+-isotherm. 

Remark.  The  hypotheses  of  Lemma  3.1  have  a  simple  interpretation.  For 
example  in  (i),  (I)  says  the  signed  area  between  the  chord  joining  (w_, 
p(w  ,0  ))  and  (w+(O,0  ),  p(w+(0,8  ),  8  ))  and  the  graph  of  p(w  ,0  ) 
between  w_  and  w+(0,8  )  is  negative.  U(0  )  is  that  positive  value  of 
U  so  that  the  signed  area  between  the  chord  joining  (w_,  p(w  ,8  ))  and  the 
graph  of  p(w,0  )  between  w_  and  w+(U(8  ),0  )  is  zero  (see  Figure  2).  An 
analogous  interpretation  holds  for  (ii). 

Since  the  arguments  used  in  describing  the  two  cases  w_  satisfies  (I)  and 
w+  satisfies  (II)  are  analogous  we  shall  assume  in  what  follows  that  w_  is 
given  in  the  a-phase  of  the  0  isotherm  and  (I)  is  satisfied. 

Lemma  3.1  assures  the  existence  of  an  isothermal  solution  of  (2.6; 

•  * 

U  (0_),o)  connecting  (w  ,0,0  )  to  (w+(U  (8  ),O,0  ).  Unfortunately  this 

solution  does  not  satisfy  the  full  system  (2.6;U  (8  ),u )  when  u  /  0.  The 

* 

reason  is,  of  course,  that  the  quantity  in  braces  in  (2.6»U  (0  ),U;c)  would 
have  to  be  identically  zero  along  this  isothermal  solution  and  it  is  not. 


* 

We  also  know  however  that  if  u  is  small  and  U  is  near  U  (8  ) 

* 

solutions  of  (2.6;U,m)  stay  close  to  solutions  of  (2.6;U  (0  ),0)  provided 
w*  or  v'  is  not  small.  If  w'  and  v*  are  small  and  V  >  0,  then  the 
"slow"  system  (2.6;U,U;c)  becomes  "fast"  relative  to  the  "fast"  system 
(2.6;U,U;a,b).  A  singular  solution  of  (2.6)U,u)  will  consist  of  alternating 

solutions  or  solution  segments  of  the  two  systems  (2.6>U  (0  ),0)  and 

* 

( 2. 6; U  (0  ),y;c),  the  latter  being  defined  where  w*  *  v'  *  0,  i.e.  on 
w  =  g(9)/  a  solution  of 

-  U*(0_)2  (w-w_)  -  p(w,0)  +  p(w_,0_)  -  0  .  (3.1) 

Our  hope  is  that  if  we  can  construct  a  singular  solution  of  (2.6;U,U) 
connecting  the  state  w  =  w_,  v  ■  0,  0  =  0  in  the  (liquid)  a-phase  to  a 
state  w  =  w+,  v  =  0,  0  =  0+  in  the  (vapor)  B-phase  we  can  find  a  true 
solution  of  (2.6;U,U)  that  does  the  same  thing,  provided  M  is  sufficiently 
small. 

From  Lemma  3. 1  we  know  the  first  piece  of  our  singular  solution  will  be  a 
curve  in  the  0  =  0  ,  v  >  0  half  plane  connecting  (w  ,0,0  )  to 

(w+(U  (0_),0_),O,0_).  According  to  our  program  the  next  piece  of  our  singular 

* 

solution  will  be  to  follow  the  flow  given  by  (2.6;U  (0  ),U;c)  when 
w  =  g(0).  Examination  of  the  van  der  Waals  isotherms  (Figure  1)  for  differing 
values  of  0  leads  to  the  graph  of  g  shown  in  Figure  3. 


I 


W+(U  (0_),0_) 

Pigure  3 

To  follow  the  solution  of  (2.6;U  (0  ),li>c)  along  w  *  g(0)  means  we 
are  studying  the  scalar  equation 

* 

»'  -  W  ~  <ej  {-(e(g(6),0)  -  ej  -  (g(0)  -  w_)p_ 

*  ,  (3.2) 

U  (0_r 

+  — —  (g(0)  -  »_)  /  . 

As  the  first  part  of  the  singular  solution  took  us  from  (w  ,0,0  )  to 

* 

(w+{U  <0_),0J,O,0J,  it  is  necessary  to  study  behavior  of  (3.2)  only  on  the 
w  >  w^  branch  of  g(0)  where  g'(0)  >  0  (Figure  3). 

Crucial  to  the  study  of  (3.2)  is  knowledge  of  equilibrium  point 
positions.  The  equilibrium  points  of  (3.2)  are  found  at  the  intersections  of 
the  graph  w  ■  g<0)  and  w  -  1(0)  where  w  ■  4(0)  is  the  solution  of 


U*(u  )2 


-  (e(w,8)  -c)-p(w-w)+ 


(w  -  w  )  -  0  . 


(3.3) 


If  we  substitute  our  constitutive  relation  (1.8)  into  (3.3)  and  use  the 


fact  that 


*  „  -p(w  (U  (6  ) , 8  ))  +  p 

u  (8  )2  - ^ — = — - - : 

W+(U  (0J,8_)  -  W_ 


we  find  w  »  4(0)  must  satisfy 

cv(e-0_)  -  a(-~  -  -^J 
w  w 


( - r 


w+(U  (8  ) ,8  )  -  w 


-)  (p_  (w+(U*(8_),0J  - 


*  w  -  w 

P(w+(U  (8_),8_),0_)  . 


(3.4) 


(3.5) 


Differentiation  of  (3.5)  with  respect  to  w  yields 
d6  -2a 


v  dw  w3 


T— L 


w+(u  (8_),ej  -  w 


-)  (P  (w  (U  (0  ) ,  0  )  -  w) 


(3.6) 


+  p(w+(U  (®_)»0_)  (w  -  w_)l  . 

So  if  the  8  isotherm  is  as  3hown  in  Figure  2  with  p(w+(U*(0  ),0  )  >  0  we 
d0  * 

have  —  <  0  for  w_  <  w  <  w+(U  (0_), 0_)  and  w  ■  4(0)  is  a  decreasing 
function  of  8  on  the  range  w_  <  w  <  w+(u*(0  ),0  ).  From  (3.6)  we  then  see 
if  cy  is  sufficiently  large  the  curves  w  ■  g(0)  and  w  *  4(0)  will 
intersect  at  three  points  in  w  -  0  plane  as  shown  in  Figure  4  giving  the 
equilibrium  points  of  (3.2).  Thus  (3.2)  has  an  equilibrium  point  8  ■  0+  so 
that  w+  =  g(0+)  >  w^  and  g'(0+)  >  0. 


’•V 


Figure  4 


w 

Y 


w+  w+(U  (0_),9_) 


I 


As  noted  in  the  preceding  paragraph  we  need  p(w+(U*(0  ),8  )  >  o. 
Examination  of  Lemma  3.1  shows  that  this  will  be  implied  by 
Assumption  U(8  ): 

For  U(8_)  as  given  in  Lemma  3.1  we  have  p(w+(TJ(8  ),8  )  >  o. 

We  also  note  the  above  result  depends  on  cy  being  sufficiently  large. 
Otherwise  the  curve  A<0)  might  fail  to  intersect  g(0)  except  at  the  point 
For  this  reason  we  make  the  basic  constitutive  assumption 
Assumption  cys 

cy  is  large  enough  to  (i)  force  l,  g  to  intersect  in  three  points  as 
shown  in  Figure  4  and  (ii)  yield  p(w+,0+)  >  0. 

Actually  (ii)  follows  by  continuity  from  p(w+(U*(0  ),6  )  >  o  when  cy 


is  large 


iV  IV 

Lemma  3.2.  The  equilibrium  point  8+  (where  w+  »  g(0+)  >  w^)  of  (3.2)  is 
asymptotically  stable. 

Proof.  From  Taylor  *s  theorem  we  know 

* 

8*  -  u  ~  (9J  {-  VvV  -  p(w+,9+)}g'(0+)(e  -  0+) 

* 

-  v  ~  (9 J  ee(w+<8+)(8  -  0+)  (3,8) 

+  0{ |8  -  8+|2)  . 

We  note  first  that  «  cy  >  0,  e ^  >  0.  Secondly  Assumption  cy 

~  —  w 

implies  p(w+/0+)  >  0.  So  8+  is  asymptotically  stable.  In  fact  the  full 
phase  flow  of  (3.2)  is  shown  in  Figure  4  by  the  arrows  on  the  graph  of  g. 

If  we  combine  Lemmas  3.1  and  3.2  we  see  we  have  constructed  a  singular 
solution  of  ( 2. 6;U,u ).  The  singular  solution  consists  of 
(i)  an  isothermal  solution  of  (2.6;U*(8_),0)  running  from  (w  ,0,8  )  to 
(w^(u  (8_),8_),0,8^)  in  the  v  >  0,  8*8  half  plane; 

(ii)  a  solution  of  (3.2)  traveling  on  the  graph  of  g  in  the  v  «  0  plane 
from  the  point  (w^(u*(8_), 8_), o, 0_)  to  the  equilibrium  of  (3.2)  0+. 

Pictorially  the  singular  solution  is  represented  in  Figure  5. 
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Figure  5.  The  singular  solution 


A®  w_  i®  in  the  a-phase  of  the  9_  isotherm  and  w+  is  in  the 
0-phase  of  the  0+  isotherm  we  see  this  singular  solution  does  indeed  connect 
a  and  0  phases. 

Our  goal  now  is  to  show  there  is  a  true  solution  of  (2.6?U,l>)  for  U 
near  U  (9_)  and  |i  small  with  the  same  connecting  properties  as  the 
singular  solution.  Before  doing  this  we  recall  some  results  from  the  Conley- 
Easton  theory  of  isolating  blocks. 


[ 


4.  Conley-Easton  theory  of  isolating  blocks 

In  this  section  the  Conley-Easton  isolating  block  theory  [2]  is 
outlined.  The  presentation  given  here  is  the  same  as  that  given  by  Carpenter 
(3]  modulo  a  small  generalization  (Corollary  4.6). 


Consider  a  system  of  ordinary  differential  equations 

y  -  F(y) 


(4.1) 


IN  N 

where  F  6  C  (ft ;R  )  for  ft  an  open  connected  set  contained  in  R  .  For  all 

initial  data  6  ft  we  shall  assume  (4.1)  has  globally  defined  unique 

solutions  y(t,yQ),  y(0,yQ)  *  yQ  lying  in  ft. 

Definition  4.1.  B  is  a  block  for  (4.1)  if  there  exists  C1  functions 


fj.fj,,..,^  i  so  that 

Ban  f  ( [0,—) ) 
i  =  1 

is  homeomorphic  to  (0,1 JN  and  f^  S  Vf^  •  F  0  on  3b.  b+  (ths  entrance 
set)  s  (y  6  3b  s  fi(y)  *  0  and  f^y)  >  0  for  some  i).  b"  (the  exit  set) 

5  {y  6  3b  :  r(y)  =  0  and  f^y)  <  0  for  some  i). 

Example.  A  block  about  the  saddle  point  y1  »  y2  «  0  for  the  system 


*1  "  y2  ' 


y2  -  y,  • 


(4.2) 


mt  b  *  1 1 y  1 1  +  | y 2 1  <  i).  b  is  a  block  with  f1(y1,y2)  =  1  -  y1  -  y2, 
f2(y1'y2)  “  1  +  Vi  +  Y2*  etc*  ^l(yry2)  “  "yl“y2  *  -1  on  fi(yi'y2)  =  °' 


so  b 


±  f“  (o)  n  b.  fjty^yj)  =  y,  -  y2  *  i  if  f2(yi,y2)  3  °» 


b+ f ^  ( 0)  n  B.  An  illustration  is  given  in  Figure  6. 


'  U 


Figure  6 


Definition  4.2.  If  f  is  a  C1  function  from  Q  to  E  and  y 
y(t,yQ)  crosses  {y  >  f(y)  -  c)  transversely  if  T  =  inf {t  >  0; 
f(y(t,yQ))  ■  c>  <  •  and  f(y(t,yQ))  /  0. 

T*(y0)  (the  time  needed  for  yQ  to  reach  b*)  = 


sup it  >  0  s  y(t,yQ)  n  b*  is  empty}  if  yQ  ^  b* 


♦*(y0>  5  y<**»y0>  if  0  *  ?*%>  <  "  • 

♦*{y  >  is  the  first  point  of  y(t,yQ)  in  b*  . 


D*  contains  the  set  on  which  is  defined  and  continuous*  If  y  C  D*, 

y(t,yQ)  crosses  ly  ;  f^y)  *  c}  transversely  at  ^(Yg)  for  some  i. 

In  the  example  given  above 

d+  «  (y1,y2  ;  y*  -  y2  <  ’>  n 

ty^Yj  ;  Y1  <  0  *nd  y2  >  °*  or  Y1  >  0  and 

y2  <  0>  -  B  , 

d"  =  (D+  u  B)  -  b"  U  S(< 0,0>)  . 

Here  S(<0,0>)  is  the  stable  manifold  of  (4.2)  running  into  (0,0). 

We  now  state  the  following  results.  Proofs  may  be  found  in  [3],  [14]. 
Lemma  4.3.  (Continuity  of  maps  defined  by  a  block).  If  B  is  a  block, 

T*  and  are  continuous  on  D*. 

Consider  the  parametrized  system 

y  -  F(y ,o)  (4.3)0) 

where  F  6  C1 (12  *  £  ;  RN)  for  E,  <2  open  connected  sets  contained  in  RN, 
1^  respectively. 

*  * 

Lemma  4.4.  If  B  is  a  block  for  (4.3»o  ),  0  6  E,  there  exists  a 

*  *  t 

neighborhood  E  of  o  ,  E  c  E,  so  that  B  is  a  block  for  (4.3 ;0)  when 
o  6  E*. 

Theorem  4.5.  Assume  there  exists  a  block  B  for  (4.3)0)  for  all  o  6  E 

with  the  following  properties. 

(i)  There  is  a  equilibrium  point  y1  of  (4.3)0),  y^  ^  B. 

(ii)  There  is  a  path  {<zg,as>  :  0  <  s  <  l}  c.  D+  *  E  so  that  zg  6  U< y1 > 

(the  unstable  manifold  leaving  y1 /  Zq,Zj  6  D“)  ♦  °  ♦+(zq»°0^» 

<p~  o  p  (a^;0  )  are  contained  in  distinct  components  of  b”. 

Then  for  some  o  one  positive  semiorbit  beginning  at  $+ (z  ;0  )  is 
s  s  s 

contained  in  B. 
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Since  Iheorem  4.5  is  slightly  different  from  the  result  given  in  [3) 
(Theorem  1.5  of  [3]),  we  present  a  proof  for  completeness. 

■f  it  x 

Proof  ofTheorem  4. 5.  Since  £  is  continuous  it  maps  a  path  u  in  D  to  a 
path  in  b+.  If  $+(II)  is  contained  in  D~,  $  °  4+(H)  is  a  path  in  b“. 

Thus  $+(II)  is  not  contained  in  D~  if  the  end  points  of  $+(II)  are  mapped 
by  4  to  distinct  components  of  b~.  In  this  case,  one  positive  semiorbit 
beginning  at  $+(II)  is  contained  in  B. 

Corollary  4.6.  Assume  the  hypotheses  of  Theorem  4.5  hold.  Assume  in  addition 
there  exists  a  C1  function  V  :  RN  ♦  R  with  the  property  that 
V(y)  =  grad  V(y)  •  F(y ,a)  <  0  for  y  6  B,  a  6  E  . 

Let 

S  *  {y  in  B;  V(y)  =  0} 

and  M  be  the  largest  invariant  set  in  S.  If  M  consists  only  of 

equilibrium  points  of  (4.3;c),  then  for  some  0  (4.3 ;o  )  possesses  a 

s  s 

heteroclinic  solution  connecting  y.,  to  y2  (y^-0*)  ■  y1 ,  y(+*)  *  y2) 

where  y0  is  an  equilibrium  point  of  (4.3>0  )  lying  in  B. 

‘  s 

Proof.  Prom  Theorem  4.5  we  know  the  existence  of  an  orbit  which  leaves  y1 
and  eventually  enters,  never  to  leave,  B.  The  function  V  is  a  Liapunov 
function  in  the  sense  of  LaSalle.  So  LaSalle's  invariance  principle  (Theorem 
1.3,  Chap.  X  of  [15])  implies  this  orbit  approaches  equilibrium  point  in  B 
as  t  ♦  ». 
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5.  Construction  ->f  the  connecting  orbit. 

In  this  section  we  use  theory  of  Section  4  to  prove  the  existence  of  an 
orbit  connecting  the  equilibrium  point  (w  ,0,9  )  of  (2.6;U, w)  to  another 
equilibrium  point  (w+, 0,9+)  where  w-  is  in  the  o-phase  of  the  9 
isotherm  and  w+  is  in  the  9-phase  of  9+  isotherm.  We  shall  need  the 
following  lemmas. 

Lemma  5.1.  For  9  =«  9  the  equilibrium  point  (w+(U  (9  ),9  ),0)  of 

W*  a  V 

*  .  2  <5*U 

Av'  =  -0  (8  )v  +  0  (9  )  (w  -  w  )  +  p(w, 9 )  +  p(w  ,9  ) 

is  a  saddle.  Furthermore  for  9  near  8  the  equilibrium  points  (g(8),0) 
of  (5.1)  near  (w+(U  (8  ) , 0  ),0)  are  also  saddles  and  possess  phase  portraits 


as  shown  in  Figure  7 


Proof.  Linearization  of  (5.1)  when  0  ■  6  about  (w+(U  (9_)*9_)»0)  shows 

the  eigenvalues  of  the  linearized  system  satisfy 

AX2  +  U*(0  )A  +  (U*<8  )2  +  p  (w  <U*(8  ),0  ),0  ))  =  0  . 

-  -  w  +  -  -  - 

Since  U*(0  )2  +  pw(w+(CJ*(0  ),0_),0  )  <  0  (see  Figure  2)  the  equilibrium 
point  is  a  saddle.  The  rest  of  the  lemma  follows  from  continuity  with  respect 
to  the  parameter  9. 

Lemma  5.2.  For  9  near  0  (5.1)  possesses  a  block  BQ  as  shown  in  Figure 

7  where  the  lengths  of  the  sides  of  Bg  are  independent  of  Bg. 

Proof.  Lemma  1 . 2  of  13]. 

Lemma  5.3.  Define  B  =  U  BQ 

0  -v<0<0 
+ 

where  0+  -  v  is  shown  in  Figure  5.  Then  B  is  a  block  for 
(2.6;  U*(9_),u). 

Proof.  By  Lemma  3.2  solutions  of  ( 2.6;U* (0_), U)  are  such  that  the  top  and 

bottom  of  B  are  entrance  points.  We  know  the  sides  are  points  of  exit  or 

entrance  from  Lemma  5.1. 

Lemma  5.4.  For  U  near  U  (0  ),  B  is  a  block  for  (2.6;U,li). 

Proof.  Lemma  4.4  applies. 

Lemma  5.5.  For  P  sufficiently  small  there  exist  Uy  >  0  and  a  solution  of 
( 2.6;U^,U)  with  w(-»)  »  w  ,  v(-»)  «  0,  8(-*)  =  9_  so  that  (w,v,9)  is  in 
B  for  all  C  sufficiently  large.  Furthermore  Uy  ♦  U  (0_)  as  M  ♦  0+. 
Proof.  When  V  *  0  linearization  of  (2.6;U  (9_),0)  about  (w_,0,9_)  shows 
(2.6;U*(0  ) ,0)  has  a  one  dimensional  unstable  manifold  exiting  from 
(w  ,0,9  )  (see  Figure  8).  By  continuity  with  respect  to  parameters,  if  U 
is  near  U*(9  )  and  9  is  sufficiently  small,  there  is  a  one-dimensional 
unstable  manifold  U(U,u)  for  (2.6;U,M)  exiting  from  (w  ,0,9  ).  Consider 
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an  arc  A*  lying  in  (j(U,y)  (A*  will  lie  near  the  arc  A  shown  in  Figure 

9).  As  U  runs  £rom  values  slightly  greater  than  U  (0  )  to  values  slightly 
* 

less  than  U  (6  )  we  see  by  continuity  with  respect  to  the  parameters  U, u 
that  the  characteristic  features  of  Figure  9  are  preserved  for  the  full  system 
(2.6;  U,M)  if  M  is  small.  Namely  the  end  points  of  A*  will  be  carried 
into  distinct  components  of  b_  6  B.  Since  Lemma  5.4  tells  us  B  is  a  block 
for  (2.6;U,ii)  for  all  U  near  U  (0_),  we  know  by  Theorem  4.5  that  for 
each  u  sufficiently  small  there  exists  so  that  (2.6;U^,y)  possesses  a 

solution  exiting  from  (w  ,0,6  )  which  enters  and  remains  in  B.  The  proof 

*  „ 

implies  may  be  chosen  such  that  ♦  U  (0_)  as  y  ♦  0+. 

Lemma  5.6.  (i)  Along  smooth  solutions  of  (2.2)  the  following  equality  is 

satisfied 

*\  *  ux<mx  -  *“2“xx>  *  “6xx  •  <5-2> 

(ii)  Along  solutions  of  (2.6|U,|i)  the  following  equality  is  satisfied 
-  6n'  *  0W'2  +  Aw*w"  +  ~  S"  .  (5.3) 

**  a  a 

a  3\b 

Proof,  (i)  Consider  the  specific  free  energy  ♦  with  p  =  -  n  -  -  yg. 
Then  since  <|>  *  e  -  0n  we  have  to  ■  e  -  0  n  -  6n  and  hence 

v  v  W  L 

A  A 

♦  w  +  <|>ft0  »  e  -  0  n  -  0n. .  This  implies  e  +  pw  »  0n  .  From  (2.2;b)  we 

W  t  o  t  t  t  t  t  t  t 

find  et  +  pu^  =  0nt»  Now  use  (2. 2?c)  to  complete  the  proof. 

(ii)  Substitute  u  =  u(C),  w  ■  w(C),  0  »  0(C),  n(C)  *  n(w(c),0(C))  into 

(5.2). 
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Figure  8 

Theorem  5.7.  (i)  (Compression  Wave)  Assume  w  ,0  ,u  are  given,  w_ 

satisfies  (I)  of  Lemma  3.1,  and  the  constitutive  Assumptions  cy  and  U(6  ) 
are  satisfied.  Then  for  u  sufficiently  small  there  exists  >  0  and  a 
solution  of  (2.6;UJ1,U)  with  w( -*)  ■*  w  ,  v( -•)  «  0,  0  ( -“)  =  0  ,  w(+<»)  » 
w+(UiJ),  v(+*°)  »  0,  0(-t")  ■  0+(U^)  <  w+(U^)  is  *,n  t*ie  ®“P^ase  ®+(U^) 

isotherm.  Furthermore  ♦  0  (0_)  >  0,  w+(U^)  +  w+»  0+(Uy)  +  0+  <  0  as 
U  ♦  0+. 

(ii)  (Expansion  wave)  Assume  w+,0+,  u+  are  given,  w+  satisfies  (II) 
of  Lemma  3.1,  the  constitutive  Assumption  cy  is  satisfied  and 
p(w+,0+)  >  0.  Then  for  V  sufficiently  small  there  exists 
U  <  0  and  a  solution  of  (2.6;U  ,li)  with  w(-*»)  ■  W*(U  ),  v(-®»)  ■  0, 


0<-*»)  -  8*0^)  >  0+,  w(+*»)  -  w+,  v(+“)  -  0,  0(+»)  -  0+>  w*(Uy)  is  in 
the  a-phase  of  the  0*(U^)  isotherm.  Furthermore  ♦  U  (0+)  <  0, 
w*(U  )  -*•  w  ,  0*(U  )  ♦  0  >0  as  u  +  0+  where  w  ,  0  are  obtained  in  a 

y  -r  -  y  .  4-  —  — 

<v  <%* 

manner  analogous  to  that  given  for  the  construction  of  w+,  0+  . 

Proof,  (i)  Define 


V(w,v,0)  =  n(w,0)  +  t-  (e(w,0)  -  e_)  -  p_(w  -  w_)  +  y-  (w  -  w_)2} 


From  (5.3)  it  follows  that  along  solutions  of  (2.6;U,)i)  we  have 


•  *  -  -g-  w,2(1  +  {-  (e(w,0)  -  e_)  -  p_(w  -  w_)  +  (w  -  wj2}) 


UP 

a02 


{-  (e(w,0)  -  e_)  -  p_(w  -  w_)  +  (w  -  wj2}2  . 


Now  if  (vfw, 0)  is  in  B 

u2  2 

{-  (e(w,8)  -  e_)  -  p_(w  -  w_)  +  —  (w  -  w_)  } 
can  be  bounded  independently  of  M.  Hence  for  V  sufficiently  small  we  have 
V'  <  -const. (w'  +  0'  ),  for  some  positive  constant,  when  (v,w,0)  is  in 

B.  Hence  the  hypotheses  of  Corollary  4.6  is  satisfied  and  the  result  follows. 

(ii)  This  follows  by  proving  a  sequence  of  lemmas  analogous  to  Lemmas 
5.1  -  5.5  for  the  case  where  (w_,0_)  is  replaced  (w+,0+)  in  (2.6jU,u). 

The  same  Liapunov  function  given  in  (i)  will  work  as  along  as  we  replace 
(e  ,w  )  by  (e+,w+).  A  pictorial  representation  of  the  singular  solution  for 


this  case  is  given  in  Figure  9. 
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/b  w  (u  (0  ),e  )  w 

' 

f 

Figure  9 

Remark.  While  the  special  case  of  parameters  (2.1)  greatly  simplifies 
calculations  results  similar  to  those  given  here  are  valid  for  the  general 
Korteweg  relation  (1.5;a).  Specifically  consider 

U(w,0)  =  6p0(w,8)  , 

C(w,0)  =  62Cq<w,8)  , 

D(w,8)  -  «2Dq(w,8)  ,  (5.4) 

h  -«8X  . 

b  =  q  =  0  . 

where  5  is  a  small  positive  constant,  a-1  -0(1)  in  6,  mq(w, 8)  >  "y  >  o 
Cq(w,6)  >  C  >  0,  u,  C  positive  constants,  and  CQ,  DQ  are  smooth 

functions.  In  this  case  (2.6;U,U)  is  replaced  by 


w' 


V 


C0(w,8)v*  -  -UM0(w,fl)v  (5.5;U,6) 

-  U2(w  -  W  )  -  p<W, 0)  ♦  P_  +  Dq(W,  0)v2  , 

0*  ■  i.  u  {-  (e(w, 0)  -  e_)  -  p_(w  -  w_) 

+  -j  (w  -  w_)2}  , 

where 

r  X  -  Ut 

5  -  — j -  . 

Instead  of  Theorem  5.7  we  can  obtain  the  following  more  general  result. 
Theorem  5.8.  (i)  Assume  w_  is  given  in  the  a-phase  of  the  0  non¬ 

monotone  isotherm  and 

w  <0,0  )  w  <0,0  ) 

<!')  Jw  '  exp  (2  fx*  "  Do<C,0_)Co(C,e_)"1dC)Co<X,8_)_1f_(X|O)dl  <  0  , 

and  there  exists  U(0  )  so  that 

w  <U<0  ),0  )  w  (0,0  )  ,  , 

Jw  "  exp<2  fx  ~  Do<C,0_)Co<e,0_)”,dC)Co(X,0_)"1f_a,U)dX  -  0  . 

Also  assume  constitutive  Assumptions  cy  and  U(0_)  are  satisfied.  Then  for 
6  sufficiently  small  the  conclusion  of  Theorem  5.7  (i)  holds  for 
<5.5, 0,6). 

(ii)  Assume  w+  is  given  in  the  0-phase  of  the  0+  non-monotone 

isotherm  and 

w  w  (0,8  ) 

UI')  /w  (0,e  )exp(2  V  D0K,e+)C0{5,V"  «)Co(X,0+)_1f  +  (X,O)dX  >  0 
and  there  exists  U(0  +  )  so  that 
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W  w  (0,0  ) 

/  exp  (2  /  "  Dn(C,ei)Cn(e,8  )"'dt)Cn(X,e  )■  f  (X;U)dX  -  0  . 

wJU<0+),8+)  A  0+0  +  0  +  + 

Also  assume  conti tutive  Assumption  cv  is  satisfied  and  p(w+,9+)  >  0.  Then 
for  6  sufficiently  small  the  conclusion  of  Theorem  5.7  (ii)  holds  for 
(5.5;  0,6). 

Sketch  of  proof.  The  main  idea  is  to  use  the  generalization  of  Lemma  3. 1 
given  by  Lemma  3.8  of  [1]  to  do  the  isothermal  connections.  After  that 
everything  can  be  done  as  before  except  that  we  now  use  the  Liapunov  function 
V(w,v,8 )  =*  n(w,8)  +  ^  {-  (e(w,8)  -  e_)  -  p_(w  -  w_)2 


02,  .2l  C0(W'9>  fv2, 

+  __  (w  -  W_>  }  +  - g -  [— J 


In  this  case  we  find  along  solutions  of  (5.5;U,6) 


V*  -  -w,2{-  UU^w^) 


CO<W'0)  6* 


3CO(W'0)  Q.  a  e-2  3co<w'6>  w.3  Do(w'6)  3 

- 575 -  ~  T  —  +  - -  TTS-  +  - H -  W 


W  ~ZZ 2 


5w  20 


As  before  the  0*  terms  in  the  expression  in  parentheses  can  be  made  small 
for  6  sufficiently  small.  The  w'3  terms  also  cause  no  difficulty  as  they 
are  dominated  by  w’*  for  B  with  sufficiently  small  sides.  Hence  we  again 
find  V'  <  -const. (w' 2  +  8’2)  when  (v,w, 0)  is  in  B. 

Remark  5.9.  We  note  that  as  M  +  0+  the  traveling  wave  solution  given  by 
Theorem  5.7  (i)  approaches  the  functions 


w(X,t) 


w 

+ 

,  u(x,t) 
w 


u  — U  (0  ) (w^  -  w  ) 
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X  >  U  (0  )t 


X  <  U  (9  )t  . 


(5.6) 


(Here  u(X,t)  is  determined  by  (2.5;b).) 

A  similar  statement,  of  course  holds  for  Theorem  5.7  (ii).  The  point, 
though  is  that  (5.6)  is  not  a  distributional  solution  of  the  inviacid  system 


(5.7) 


ufc  ■  -p(w,0)x  , 


Wt“UX  ' 


Et  -  “(up)x  +  °®Xx  * 


Mathematically  the  reason  is  that  (5.7)  cannot  possess  distributional 
solutions  with  a  jump  in  9?  u,w  piecewise  constant  and  (5.7jc)  implies  6^ 
is  a  locally  finite  Borel  measure  and  hence  9  is  continuous.  This  has  been 
noted  by  Dafermos  [16];  a  physical  explanation  is  given  in  {[13),  Vol.  II,  p. 


Hence  for  the  problem  considered  here  the  internal  structure  is 
crucial.  Solutions  of  the  viscous  problem  do  not  converge  to  solutions  of  the 
inviscid  problem  as  u  ♦  0+.  This  is  a  reflection  of  the  h  »  a®xx  and 
<**■0(1)  in  U  assumptions  and  not  the  unusual  non-hyperbolic  nature  of 
(5.7). 

Remarks. 10.  We  note  that  Theorem  5.7  (i)  says  that  the  states  in  a-phase 
determines  both  an  upstream  state  in  the  0-phase  and  a  speed  of  propagation 
>  0  so  that  a-phase  (liquid)  travels  into  the  0-phase  (vapor).  This 
stands  in  contrast  to  classical  fluid  dynamics  (pw  <  0)  where  both  the  state 
on  one  side  of  a  shock  and  the  speed  of  propagation  determine  the  state  on  the 
other  side  of  the  shock.  A  similar  interpretation  can  be  given  for  Theorem 


5.7  (ii) 


Remark  5.11.  It  may  be  valuable  to  note  that  for  u  small  Theorem  5.7 

*  * 

predicts  that  the  wave  speeds  U  (0  )  in  (i)  or  U  (6+)  in  (ii)  will  provide 
a  good  approximation  to  the  true  wave  speeds.  Thus  while  the  isothermal 
equation  (2.6 ;(J,0)  yields  the  wrong  equilibrium  states  to  which  a  transition 
is  made  it  does  yield  a  good  approximation  to  the  correct  speed  of  transition 
(Lemma  3.1). 

Remark  5.12.  While  it  is  not  the  purpose  of  this  paper  to  predict  any 
particular  set  of  experimental  data  it  seems  worthwhile  to  note  the 
observations  of  Dettleff,  Thompson,  Meier  and  Spectman  [17].  In  that  paper 
they  noted  the  ability  to  produce  a  wave  which  yields  complete  liquefaction  of 
a  superheated  vapor,  i.e.  a  complete  transition  from  gas  to  liquid.  The 
liquids  used  in  their  study  were  of  "retrograde  type"  (their  terminology)  in 
that  they  possessed  high  specific  heat  at  constant  volume  (e.g. 
fluorocarbons).  They  also  noted  an  increase  in  temperature  from  vapor  to 
liquid  phases.  Remarkably  all  these  observations  are  consistent  with  the 
theory  presented  in  this  paper. 

The  results  of  [17]  are  for  a  compression  wave  only.  I  know  of  no 
experimental  results  showing  the  liquid  to  vapor  expansion  wave  predicted  by 
Theorem  5.7  (ii).  However  within  the  range  of  solids  (e.g.  iron)  exhibiting 
phase  transitions  both  expansion  and  compression  waves  have  been  observed  [13; 
Vol.  II,  p.  751].  The  compression  waves  produced  an  increase  in  temperature, 
the  expansion  waves  a  decrease.  Again  this  is  consistent  with  the  conclusions 
of  Theorem  5.7. 
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